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1 Introduction 



In the Euclidean version of field theory, we are interested in computing the Schwinger functions 

of a theory. In order to obtain these functions, Parisi and Wu introduced the stochastic quanti- 
zation [1]. This formalism was introduced as an alternative quantization scheme, different from 
the usual canonical and the path integral field quantization, based in the Hamiltonian and the 
Lagrangian respectively The method starts from a classical equation of motion, but not from 
Hamiltonian or Lagrangian, and consequently can be used to quantize dynamical systems with- 
out canonical formalism and therefore it is useful in situations where the others methods lead to 
difficult problems. 

The main idea of the stochastic quantization is that a d-dimensional quantum system is equiv- 
alent to a {d + l)-dimensional classical system which undergoes random fluctuations. Some of 
the most important papers in the subject can be found in Ref. [2]. A brief introduction to the 
stochastic quantization can be found in the Ref. [3] and Ref. [4]. See also the Ref. [5]. 

In a previous paper [6], we studied the stochastic quantization of a self-interacting scalar field 
theory, assuming a non-Markovian process, modifying the Langevin equation by introducing a 
memory kernel [7] [8] [9]. We have shown that although a system with a stationary, Gaussian, 
non-Markovian Langevin equation with a memory kernel and a colored noise converges in the 
asymptotic limit of the Markov parameter r to the equilibrium, we obtain a non-regularized 
theory. 

In this paper we would like to continue to investigate the virtues of this non-Markovian stochas- 
tic quantization method, now employed in the case of a topological field theory. One of the peculiar 
features within this kind of theory is the appearance of a factor of i in front of the topological 
action in Euclidean space. Since the topological theory does not depend on the metric of space- 
time, the path integral measure weighing remains to be e*'^ even after the Wick rotation. Another 
feature of a topological action is that it is the integral of a density which is not bounded from 
below in Euclidean space. So, if one attempts to use a Markovian Langevin equation with a white 
noise to quantize this theory, one will find serious problems if the factor of i is ignored. This 
Langevin equation will not tend to any equilibrium in the large r limit. So, in this sense, the use 
of a Langevin equation with a complex action [10] becomes essential for stochastically quantizing 
a topological action [11] [12] [13]. 

There is, in the literature, an approach to solve the above mentioned convergence problem. 
Studying the purely topological Chern-Simons theory, Ferrari etal, introduced a non-trivial kernel 
in the Langevin equation [14]. On the other way, Wu et al [15] showed that the Langevin equation 
for a Maxwell-Chern-Simons theory converges to the usual equilibrium result without the need to 
introduce such kernel. Their method, however, only works in the case where the Chern-Simons 
coefficient is real. 
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We show in this paper that, if one uses a non-Mar kovian Langevin equation with a colored 
random noise, tliis convergence problem may be solved in a different way. We will apply this 
approach to three-dimensional abelian Chern-Simons theory and prove that we obtain convergence 
towards equilibrium even with an imaginary Chern-Simons coefficient. To simplify the calculations 
we assume the units to be such that h — c— 1. 



2 Stochastic quantization of abelian Chern-Simons theory 

Let us consider the following action for the three-dimensional Maxwell- Chern-Simons theory, 
in Euchdean space: 

Jd^x (^-^A^{x){-AS^, + d^d,)A,{x) - i-^e^,pA^{x)d,Ap{x)^ , (1) 

where A is the three-dimensional Laplace operator. At the end of our calculations we set e — cxd to 
obtain the results for the purely topological theory, as discussed in Ref. [15]. Notice the factor of 
i in front of the topological term, as mentioned before. In order to obtain the Schwinger functions 
o the theory, let us use the stochastic quantization method. Let us introduce a non-Markovian 
Langevin equation given by 

d r SS 

—A^{t,x)^-J^ dsMA{T - s)jj-^\A^{x)=A^is,x) + Vn{r,x), (2) 

where M\{t — s) is a memory kernel and A is a arbitrary parameter. We will have, from Eqs.(l) 
and (2), in momentum space: 

^A,{T,k) = -'^(^6,^-'^y^^dsMA{r-s)A^{s,k) + 

-^^ixvpkp ds Ma{t - s)A^{s, k) + r]^{T, k) (3) 

where the stochastic field rj^ij, k) satisfies the modified Einstein relations: 

(r/^(T,^)), = (4) 

and also 

( r7^(r, k) rj^ir', k') ), = 25^Mk{\t - t'\) 5\k + k'). (5) 
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For the initial condition A^It, k)\r=o = 0, it is easy to see that the solution of the Eq.(3) is given 
by: 

poo 

A^{T,k)= dT'G,,{k;T-r')r]u{T',k), (6) 
Jo 

where we introduced the retarded Green function G^,^(A;,r), which satisfies: 



K /"T 

--^^^^P<y^'y ds Ma{t - s) Gp^{k, s) + S,^„ 5{t), (7) 
47r Jo 

for r > and G^i.(A;, r) = for r < 0. 

To proceed the calculations, let us introduce the Laplace transform of the Eq.(7): 

zG,.{k,z) = -^(^S,,-^^M4z)G,.{k,z) + 

-^^hp<tKMk{z) Gp^{k, z) + 5^^, (8) 



where: 



Mj,{z)^ / drMA{T)e-'\ (9) 
Jo 

For the result without memory (or, formally, when M\(t) — > 5{t)), we have, from Eq.(7): 



whose Laplace transform reads: 



-^(^^pakaGp^{k, t) + 5p,^5{T), (10) 



^2 / h h 



zGf,^{k,z) = -^[^,^p — -grjGpA^^^) 



— -^^lj,pakaGpi,{k, z) + 5fj,i,. (11) 

Note the similarity between Eqs.(8) and (11). The solution to Eq.(lO) is given by [15] 

kj^ki, . j f IT k^ky 



■.,„(^,.))exp(-^.). (12) 



^^""^ k 
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whose Laplace transform is: 



^) = + ^-^ — :2 . .2 ^ (13) 



+ +(£) 



j ' V 47r 



Comparing Eqs.(8) and (11), it is trivial to obtain the analog of Eq.(13) with memory: 

^ _ UjjmA k ( £L 



where: 



and 



1 _ M^{z) 



~I2 g-2 



(15) 



k' = acMa(^). (16) 
In the appendix A, we derive in detail the inverse Laplace transform of Eq.(14). It is given by: 

G^AKr) = (^^ + g,,G,{k,r)+g,,G2{k,r)y{T), (17) 

where the quantities Gi{k,T), i = 1,2, g^i, and g^i, are defined in the appendix A. We see that 
our Gn^{k,T) does not approach zero as r — > oo. The reason of such behavior is the presence of 
the longitudinal term which is common in the stochastic quantization of all gauge theories 
without gauge fixing and can be eliminated by a suitable stochastic gauge fixing. In spite of this, 
the presence of this term will not give any contribution to gauge invariant quantities. 

After this discussion, we are able to present the two-point correlation function. We have that 
D^^{k;T,T') is given by 

D^,{k;T,T') = {A,{T,k)A,{T',k')),= 

= S\k + k') ds ds' G^^ik, T-s) Gx.{k, t' - s') { rj,{s, k) rtx{s', k') ), 

/•OO roo 

^26\k + k') ds ds'Gf,x{k,T-s)Gx,{k,T' -s')Ma{\s-s' \). (18) 
Jo Jo 
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So, inserting Eq.(17) in the above equation and splitting the result in five different contributions 
yields: 

D^,{k; T, r') = 26\k + k') f Ji + J2 + J3 + J4 + 



where: 



Ji^ I dsj ds'^Mj,{\s-s'\), 



^0 



s')Ma{\s- 




s')Ma{\s- 


s'\), 


s')Ma{\s- 


s'\), 



and finally 



J5 = £ds£ ds'~g^x9x.G2{k;T-s)G^{k-T' -s')M^{\s-s' |). 



(19) 

(20) 
(21) 
(22) 
(23) 

(24) 



We can solve these equations by ordering the fictitious times s and s', s > s' for instance, and 
solving the integrals in s {s') in the interval [0, t] ([0, s\). We obtain for Ji, in the hmit r — > 00, 



J 1 k^ky 



A2 



(25) 



The integrals J2 and J3 can be solved by analogy with the scalar case [6]. Making the following 
replacements: 

a . .0 {a^ - yl) 



2 , 2^ « .1 , A4^ a2„. ("^ -1/?) 



ik' + m%^-{l + A'')-A%- 



(26) 
(27) 



where the subscript 1 (2) stands for the Gi {G2) case (see the appendix A), we will have, in the 
asymptotic limit r — > 00 that 



J2 



^,(1-A^)+AV+ ^2 



A2 ^A^ ^ (a_+ 7)^ 



A^ ((x + 0')^' A 



((77)2 V 4 
2 



47r 



47r 



(28) 
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and 



J. 



(a7)2 V 4 



An 



(a - 7) 



2\ 2 



(29) 



The remaining integrals J4 and J5 can be solved without any further complications [16]. Again, 
in the asymptotic limit r ^ 00, we obtain: 



Ja + J5 



/(A, ^,7) 



A^ 



(20-7) V47r 



A^ 



+ 



a + 7)2\/A4 (a -7) 



+ 



A', 2 



K 

An 



,(30) 



where: 

/(A, ^,7) 



153A^^ + A^'^ 

+^((^ + 7)^ + (^-7)^ 



18((7 + 7)^ + 17((7 - 7)^ + 9(77 

9 
2 



+ A^ 



(<7 + 7)' + (a-7)' + 



cr7^((T + 7)^ — (77(1 — 2(77) 



+ 



+AV7 



- (77(a + 7)^ -^(a + 7)^(^-7)^' 



and 



^(A, a, 7) = (9A^ + ((7 - 7)') (qA^ + ((7 + 7)') (A^ + a') (A^ + f ) . 



(31) 



(32) 



As mentioned before, the linearly divergent longitudinal term, found in Eq.(25), can be eliminated 
by a stochastic gauge fixing. Now, taking the limit e — > 00, it is easy to see that the contribu- 
tion J2 + J3 vanishes identically. Then, finally, we obtain, for the purely topological two-point 
correlation function: 



D^,{k-T,T') = 26\k + k') 



2 V A2y' + 5'(A,(7,7)V 2Q{y[) Un J ^ 



+ 



I] s, 



'p,v 



k^ky 



(33) 
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where P' = p\e^oo = x^^(i^)' v'l = ViU^oc and: 

1 Q 

/'(A, a, 7) = /le^oo = 120A^^ + -fA''Qiy[) + 9A'Q'iy[) + 

g'{A,an) = gle^oo = QAA'Q'iy[) + 32A^Q'{y[)+AQ'{y[), (35) 
Q{y[) ^ iA'y[ - {y[rj . (36) 

We see that in our last expression for the propagator remained a term proportional to the Maxwell 

transversal propagator. This is a anomalous situation, since the Maxwell contribution is absent 
in the usual purely topological Chern-Simons theory. The origin of this anomalous situation is 
the use of a non-Markovian Langevin equation. To circumvent this problem and recover the usual 
result, we have to make the following choice: 



which lead us to: 



So, if we choose: 



P' = Q'{y[): (37) 

, A (A8-4/9')'/' . . 

y[-^± 2 • (38) 

A , (A8-4/3')V2 c 

where C is a real constant and n is an arbitrarily large integer number. Inserting this latter 
equation in Eq.(A.23), we will get a cubic equation in C. Prom the usual Galois theory of radical 
solutions for polynomials [17] [18], we can always choose a real root from the three possible ones. 
So, in other words, we can always choose a real constant such that the two-point correlation func- 
tion converges to a "purely topological" term, with some minor differences from the usual one. We 
notice as well that our approach still works when k, is purely imaginary (which is mathematically 
analogous to writing = A'^ + iA'^, where A'^ is real, and taking the real part of the Langevin 
equation (3) in coordinate space). 



3 Conclusions 

In this paper we discussed the stochastic quantization for Maxwell Chern-Simons theory using a 
non-Markovian Langevin equation and examined the field theory that appears in the asymptotic 
limit of this non-Markovian process. 
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This paper is the second one of a program where it is investigated the possibihty that the 
Parisi-Wu quantization method can be extended assuming a Langevin equation with a memory 
kernel with the modified Einstein relations. To make sure that this modification can be used, one 
must first check that the system evolves to the equilibrium in the asymptotic limit. Second we 
have to show that converges to the correct equilibrium distribution. We proved that although the 
system evolves to equilibrium, in the propagator remained a term proportional to the Maxwell 
transversal propagator. This is a anomalous situation, since the Maxwell contribution is absent 
in the usual purely topological Chern-Simons theory. To circumvent this problem and recover the 
usual result, we have imposed a constraint in the parameters of our theory. 
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A Appendix 

In this appendix, we derive the retarded Green function for the diffusion problem Gij,y{k,T) 
Expanding the denominator in Eq.(A.l), given by 

e _ fcpfe^A h ( ^ 



we have: 



where: 



and: 



-+^) +(£) 



GAKz) = ^\ + - ^-^)h{z) - e,.^K{£)h{z), (A.2) 

^i(^) = (A-3) 

Piz) ^z' + 2^M4z)z + ^Miiz) + (|-)'a;^M1(z). (A.5) 
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Using the following exponential representation for the memory kernel M\{t): 

Ml{r) = ^A^(A^ I r |)"exp(-A^ | r |), (A.6) 
where A is a parameter, we will have, for the case n = 0: 

'^^'^ - — W) — ' ^^-'^ 

/.W = (A.8) 



and: 
where: 

and 



q;=^, (A.10) 

In order to get the inverse Laplace transform of Eq.(A.2), we must seek for the solutions of 
the quartic equation ^2(2;) = 0. As it is well known, a general quartic equation is a fourth-order 
polynomial equation of the form: 

2;^ + as^^ + 02-2^ + aiz + ao = 0. (A. 12) 

Using the familiar algebraic technique developed by Ferrari and Cardano [19], it is easy to show 
that the roots of Eq.(A.12) are given by: 





1 1^1^ 

-r^+2^+2^' 


(A.13) 


Z2 = 


1 1 1 

-.a, + -R--D, 


(A.14) 


2^3 = 


1 1 1 ^ 

— as R+ -E, 

4 2 2 ' 


(A.15) 


2:4 = 


1 1 1 ^ 

-^a,--R--E^ 


(A.16) 
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where: 



\ 1/2 

F{R) +Gj for i? 7^ 

N 1/2 

F(0) + H) for i? = 0, 

\ 1/2 

F{R) -G) for i? ^ 

/ \ 1/2 

iF{0)-H] fori?=0, 



F(i?) = -al - i?2 - 2a2, 

G = ^(40302 - 8ai - al)R~^, 
and Hi is a real root of the following cubic equation: 

- a2y^ + (0103 - 4ao)y + (4a2ao -a\- OgOo) = 
Therefore, the inverse Laplace transform of Ii{z) and l2{z) reads: 

zl + 2K^zl + A^^i 



{Zi - Z2){Zi - 23) (Zi - Zi) 

zl + 2A22| + A^^a 



(^2 - 2^1) (2^2 - Z^){Z2 - Z4) 
zl + 2A^zi + A^^3 

{Zs - Zi){z3 - Z2){Z3 - Z4) 
■(^4-^1) (^4 -^2) (^4 -^3)^ ' 



and 



+ ~ 



(2^1 - 2:2) (2:1 - 2:3) (^1 - Z4) 
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{Z2 - Zi){z2 - Z3){Z2 - Z4) 
A% , A* 

(2:3 - 2^1) (2^3 - Z2){Z3 - Z4) 

+7 J'' ^ t T^"^- (A-25) 

(2:4 - 2l)(24 - Z2)(2;4 - Z3) 

Now, let us study a simple convergence criterium in order that G^v{k,T) ^ as the Markov 
paramter goes to infinity, i.e., r — > 00. In this situation, the system converges to an equilibrium. 
Comparing the polynomial Vl{z) with expression Eq.(A.12), it is trivial to make the following 
identifications: ao — (3, ai — aA^, a2 — ck + and, finally, 03 = 2A^. 

For convenience, let us assume that R, defined by Eq.(A.17), does not vanish. To proceed with 
the calculations, let us introduce the following real quantities a and 7 defined respectively by 

a2-\al-y,y = {a-y,y/' (A.26) 
and 

7 ^ (a2 + 2/1 - lalY' = {a + y,- K'f'\ (A.27) 
where we used the identifications 02 = o; + A^ and 03 = 2A^. Then, we shall have: 

R = ia, (A.28) 

and 

E = i-f. (A.29) 

So, with the above identifications, it is easy to see to prove that G, defined by Eq.(A.22), vanishes 
identically. Therefore, we will have, from Eq.(A.18) and Eq.(A.19), that D = E. We also see that: 

(7^ + 72 = 2q; - A^ > 0, (A.30) 

which implies: 

^2 A2 

> (A.31) 

where we used Eq.(A.lO), which is a convergence criterium similar to the massless scalar field case 
[6]. 



11 



Thus, from Eqs.(A.13) - (A. 16), Eq.(A.28) and Eq.(A.29), we obtain the following solutions to 

n{z) = 0: 

z,^-— + -ia + -tj, (A.32) 

Z2 = -— + -la - -^7, (A.33) 

1 1 /AX 

^3 = -Y- 2^cT + -^7, (A.34) 
A' 1. 1. 

^4 = - Y - ~ -11- (A.35) 
So, from these last results, we will have, finally, for Gni,{k,T): 

G,.{k,T) = (tl^ + g^^G,{k,T) + ~g,,G2{k,T)\9{T), (A.36) 



where: 

G4k,r) . (^si„((^r) +cos((^r))e-*'. (A.38) 
and Qi^i, and appearing in Eq.(A.36) are defined by: 

9,j.u = ^u,u - V, (A.39) 

Qnu = V - n^z/, (A.40) 
with: ^ 

and 
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